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The interference between extrinsic and intrinsic losses in x-ray absorption fine structure (XAFS) is 
treated within a Green's function formalism, without explicit reference to final states. The approach 
makes use of a quasi-boson representation of excitations and perturbation theory in the interaction 
potential between electrons and quasi-bosons. These losses lead to an asymmetric broadening of 
the main quasi-particle peak plus an energy-dependent satellite in the spectral function. The x-ray 
absorption spectra (XAS) is then given by a convolution of an effective spectral function over a one- 
electron cross-section. It is shown that extrinsic and intrinsic losses tend to cancel near excitation 
thresholds, and correspondingly, the strength in the main peak increases. At high energies, the 
theory crosses-over to the sudden approximation. These results thus explain the observed weakness 
of multi-electron excitations in XAS. The approach is applied to estimate the many-body corrections 
to XAFS, beyond the usual mean-free-path, using a phasor summation over the spectral function. 
The asymmetry of the spectral function gives rise to an additional many-body phase shift in the 
XAFS formula. 

PACS numbers: 71.10.-w 78.70.Dm 



I. INTRODUCTION 

The treatment of inelastic losses in x-ray absorption 
fine structure (XAFS) has long been of interest ||, ||- 
Two types of losses are identified. Extrinsic losses oc- 
cur during the propagation of the photoelectron, and 
are caused by the creation of excitations such as plas- 
mons, electron-hole pairs, etc. Intrinsic losses refer to 
the creation of excitations by the sudden appearance of 
the core hole. The intrinsic losses are often called shake- 
up and shake-off excitations but are of the same type 
as the extrinsic, i.e. plasmons, multi-electron excitations, 
etc. Typically these excitations are observed only weakly 
in x-ray absorption spectra (XAS). The extrinsic losses 
cause a decrease in intensity in the no-loss or primary 
channel, which is usually treated phenomenologically in 
terms of a mean-free path A. Owing to the difficulty 
of quantitative calculations, these additional losses have 
usually been neglected or represented as a constant am- 
plitude factor, on the understanding that they only give 
some smooth background. 

The question of possible interference between extrin- 
sic and intrinsic losses has long been unsettled. For pho- 
toemission spectroscopy, it has been shown [Q, ||] that 
this interference is particularly important near excita- 
tion thresholds, where the losses strongly cancel. This 
cancellation results from the opposite signs of the cou- 
pling between the photoelectron and the core-hole to ex- 
cited states of the valence electrons. In particular, for 
plasmon creation at threshold, only the long-wavelength 
plasmons appear which, due to momentum conservation, 
makes the recoil of the electron in the extrinsic losses 
small. The intrinsic losses are caused by coupling to the 
core hole, which cannot recoil. As a result the extrinsic 
and intrinsic couplings become equal at threshold. The 
situation in XAFS is less clear. Fujikawa has discussed 



this cancellation using an explicit calculation of matrix 
elements and approximate closure relations ||. However, 
his discussion is limited to the vicinity of threshold, while 
we find that the cancellation effects extend over a wide 
energy region. 

In this paper we present a formal analysis of the loss 
problem, which is an extension of work by Hedin and 
Bardyszewski ||, ||, ||, together with numerical calcu- 
lations to illustrate the theory. The results here are 
formulated in terms of an effective one-particle propa- 
gator which includes both inelastic losses and interfer- 
ence effects. This propagator contains an asymmetric 
quasi-particle peak plus a broad energy dependent satel- 
lite structure. Our approach is essentially a general- 
ization of the GW approximation which, in addition to 
extrinsic losses, treats intrinsic losses and interference 
terms. The formalism also partly accounts for edge- 
singularity effects and contains corrections to the final 
state rule. Although the cumulant expansion was suc- 
cessfully used to describe intrinsic losses in valence elec- 
tron photoemission |7), we have found it more difficult 
to apply for x-ray absorption spectra. For photoemission 
(PES) the propagator for a hole and the spectral function 
A h (u) = A{lo)6{E f - lu) are needed, while in XAS we 
need the propagator for a particle and the spectral func- 
tion A p (uj) = A(lu)9(lu - E F ). The GW approximation 
for A(uj) was discussed in Ref. || (p 87, 92). The satellite 
of the hole spectral function is very strong and sharp at 
the bottom of the band, but decreases in intensity and 
broadens as one approaches the Fermi energy Ep [|[ [|. 
At Ep, the satellite for the particle spectrum is similar to 
that for the hole spectrum, and it rapidly becomes weaker 
and broader with increasing energy above Ep. The GW 
approximation places the satellite at about 1.5 lu p , rather 
than at u p away from the quasi-particle peak, as is pre- 
dicted by the cumulant approximation, and born out by 
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PES experiments 0. However, above the Fermi level 
the difference between the GW and cumulant approxi- 
mations is less pronounced. 

Our formulation clarifies the nature of inelastic losses 
in XAS and also yields semi-quantitative estimates of 
their effects, based on the electron gas approximation. In 
particular the theory yields an estimate for the reduction 
in the XAFS amplitude due to inelastic losses in terms of 
a phasor summation over the spectral function. This re- 
sults in an energy dependent reduction factor |Sq(w)| to 
the usual XAFS formula, as well as an additional many- 
body phase shift. Near the excitation threshold, we find 
that there is appreciable cancellation of strong extrin- 
sic and intrinsic losses by the interference terms. Cor- 
respondingly, the strength of the primary channel, i.e., 
the main quasi-particlc peak, increases near threshold. 
Thus the theory also explains the surprising weakness of 
multi-electron excitations in the observed XAS p0| . At 
sufficiently high energies both the extrinsic and the inter- 
ference contributions become negligible, and the theory 
crosses over to the sudden-approximation limit, where 
only intrinsic losses remain. Our theory is illustrated for 
the case of Cu metal. This system provides a good test 
case for our theory since errors in the conventional MS 
expansion and potentials are minimal and accurate ex- 
perimental data is available . 



II. THEORY 

A. Basic expressions 

The x-ray absorption coefficient can be expressed 
formally in terms of the many-electron Green's function 
G(E) = 1/(E-H + as 
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where |^o) is the iV-particle ground state of the total 
system (valence electrons and ion cores) , H the Hamilto- 
nian which fully includes electron-electron interactions, 
Eq the ground state energy (we consider for simplicity 
only the case when the temperature is zero), and ui the 
photon energy. We use atomic units m = |e| = H = 1, 
where lengths are in Bohr radii (0.529A) and energies m 
Hartrees (27.2 eV). Further 



A = E k {k\d\b)clb + he 



(2) 



is the dipole operator coupling the photon to the elec- 
tronic system, and 7 the inverse core hole lifetime. We 
have assumed that a specific core level \b) on a specific 
atom is involved. 

The standard way to proceed from Eq. (|l|) is to insert a 
complete set of interacting states \^ n ) , which recovers the 
golden rule expression when we take 7 as infinitesimal; a 



finite 7 just gives a Lorentzian broadening, i.e., 

H (w) = Y, K*n|A|*o>| 3 $(E + w - E n ). (3) 

n 

This expression contains explicit final states |$ n ) in 
which the excited electron (the photoelectron) is corre- 
lated with the valence electrons. Such states are very 
difficult to handle, so we instead take a different route 
which leads to an expression where no explicit final states 
are involved but instead, an expansion in one electron 
Green's functions, as in conventional XAS theory. Our 
formulation, however, contains shake up effects and ex- 
trinsic inelastic losses as well as interference terms, and 
not just the damping in the elastic channel. 

Some of the formal results derived here were presented 
previously in short conference reports || p| . In this paper 
we have concentrated on the shake-up effects of the core 
hole potential and on the extrinsic losses of the ejected 
electron. To have a clean picture, we have regarded the 
core electron as structureless. Thus we have not consid- 
ered interesting problems like core hole degeneracies and 
L///L/77 edges. Such problems are in principle compli- 
cated Kondo type problems. If we, however, limit the 
treatment to have an ion embedded in a solid state envi- 
ronment, and only consider the multiplct effects for the 
embedded ion , we think our treatment could straight- 
forwardly be extended to this situation. 

We take for H the standard "deep-level" Hamiltonian 
discussed e.g., by Langreth || (cf. Ref. 0, pp 645), 



H = H v +e c tfb + Vbtf . 



(4) 



Here H v describes the electrons outside the ion cores, 
i.e., the valence electrons and the photoelectron, and V 
the interactions between the outer electrons and the core 
hole. This approximation neglects virtual excitations of 
the core electrons, but takes correlations among the outer 
electrons and the potential V from the core hole fully 
into account. Neglecting core- valence exchange we have 
explicitly, 

V = -^2w(Ti), w(r)= fv(r-r')p b (r')dr', (5) 



where v(r) is the Coulomb potential, and Pb(r) the charge 
density of the core electron state "6" . The ground state 
of H l^o); is t nus a product, 



|*o) = |*o>|6>- 



(6) 



where |$q) is the correlated wave function for N v outer 
electrons and |6) is the core electron wave function, 



H v |$ ) = £ Q ° |*o) , H |* ) = E |*o) , E = £ ° + e c . 

The passive core electrons are not written out explicitly. 
The core electrons are easily eliminated, and Eqs. dll), (0) 
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and (||) give 
where 



-Im y(b\cfi\k 1 )(k 2 \d\b) 
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X ^\c kl -— 1 TT ,.- c ll^o), (7) 

E + u> - H' v + ij 2 



H' V = H V + V. 

Here and elsewhere in this paper, we will use a prime 
to denote quantities calculated in the presence of a core 
hole. The x-ray absorption is now formally given by a 
one electron expression 

M H = --Im V (b\rf\ki){ki\ gi B (w + £ c ) |fe)(fc 2 |d|6), 
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(8) 



where g ff (w) is an "effective" one electron Green's func- 
tion 



(AilSeS (w) |fc 2 ) = ($o|c fcl 



1 



-4j$o). 



(9) 



The quantity £7q is the ground state energy of H' v for iV„ 
electrons, and E c = e c + Eq — E' is the renormalized core 
electron energy. This Green's function is not of standard 
form since |<&n) is an eigenfunction of H v and we have 
H' v in the denominator. 

The theory developed so far is quite general, and can 
even account for losses, threshold singularity effects, and 
deviations from the "final state rule," i.e., the prescrip- 
tion that the XAS is given by a one-electron expression 
with dipole matrix elements between the initial core and 
final state wave functions calculated in the presence of the 
core hole. The function (u>) is always positive, and g c g 
can be written in terms of a Hermitian spectral function. 
An approximation for the "transient" Green's function 
in Eq. (^|) to describe the edge shape was given in Ref. 
|l4| , pp 674. However here we want to describe loss pro- 
cesses, and thus we have to develop different approxima- 
tions. These approximations, as discussed below, can be 
summed up in the quasi-boson representation described 
in Sec. II C. 

To get a qualitative feeling for the properties of the 
transient Green's function, we first discuss them in the 
Hartree-Fock approximation. The ground state of H v , 
|$o) then is a Slater determinant, and the intermediate 
states in Eq. (^) are Slater determinants built from or- 
bitals which are self-consistent solutions of a Schrodinger 
equation with the core-hole Hamiltonian H' v = H v + V. 
We are free to choose any complete set to represent 
the states k in the optical transition operator A = 
Ttk(k\d\b)cl.b. Here, we take the states that belong to 
H v which have the convenient property that cj. \Qq) = 



for k < kp- In H' v we single out one term hi with or- 
bitals that describe the photoelectron, and one H' v0 with 
orbitals for the rest of the system, 



H'v - h ' + H 'v0- 



(10) 



These two terms both have interactions with the core 
hole, V pc and V vc respectively, which are screened ver- 
sions of V in Eq. (^), since we consider a self-consistent 
HF solution for the core hole Hamiltonian. In terms of 
the old Hamiltonian H v = h + H v q , 



h' = h + V pc , 



H v o + V v , 



(11) 



With the h! orbitals separated from the H' v0 orbitals, 
only the one-electron operator h! can couple c\, to c k 
in Eq. (^). In such a case we can use a product space, 
c k |$o) = I'I'o) which then gives for the many-body 
XAS, 



H(u) = Im (b\dPg cS (w + E c )Pd\b), (12) 

where 



eff (w) = ( $. 



{H' v - E' ) 



(13) 



and the projection operator onto unoccupied one-particle 
states of the initial state (without a core hole) is 



(14) 



k>k f 



The separation of H' v into h! and H' v0 is clearly an approx- 
imation. However, it makes physical sense and has been 
used previously, e.g. for photoemission problems Q. In 
this expression the coupling between the photoelectron 
and the valence electrons is not present, since it is a cor- 
relation effect beyond the HF approximation. This de- 
fect is of less importance for a localized system, where 
HF theory is often quite useful. 

Let us now turn to the correlated case. A correlated 
|$o) can in principle be calculated from configuration 
interaction theory. Then |$o) is a sum of Slater deter- 
minants having differing numbers of electron-hole excita- 
tions (virtual excitations). The Slater determinants with 
virtual states close to the Fermi level have the largest 
coefficients in this expansion. For photoelectron states k 
away from this virtual cloud it is a good approximation 
to use a product space. For definitcness we will use 



|$o)|fc), k > k F 
0, k < kF, 



(15) 



being aware that this is a dangerous approximation for 
k ~ kp. As before we split H' v into one part hi which 
describes the photoelectron, and one part H' v0 for other 
excitations of the valence electrons, and we introduce core 
hole potentials as in Eq. (pd|). In addition we now also 
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have the dynamic coupling V pv between the photoelec- 
tron and the valence electrons, 



val 



( c l c h 



(16) 



The term (kil\ \ \v\ | is an antisymmetrized matrix el- 
ement of the Coulomb potential v(r), and the expectation 
value (c; c z ) is subtracted, since it is already included in 
the definition of h! . The state |$o) is an eigenfunction of 
H v q built from one-electron cigenfunctions without core 
hole potential while b! is built from eigenfunctions \k') 
with a core hole present, 



ti 



E ^W){k'\ 



(17) 



k'>k F 



The states \k') are scattering states. There is thus a 
one to one correspondence between \k') and |fc), and the 
energies are unchanged. Since h' only has terms bilinear 
in the photoelectron operators c^, and since there is a 
linear relation between the states |fc') and h! will not 
take us outside the product space |<& )|fc). 

From Eqs. ([?]) and (|l5|), we again obtain Eq. ([[2]) for 
the XAS, but with g c s(u)) in Eq. (||) replaced by 



$0 



{H' v0 -E' )-N-V vv +i 1 



(18) 



where uo n — E' n — E' a . Putting in eigenstates \k') of h', 
and taking the imaginary part, we then obtain, 

M H = J2\(^o\K)\ 2 \(k'\Pd\b)\ 2 5(uj + E c -uj n -e k ) 



k,n 

u>+E c -E F 



dw' A{l}')li {1) {lu + E c -lu'), (20) 



where the core-hole spectral function A{u>) is 
^H = El(*oK)| 2 <5(c-c„), 

and the one-electron XAS is 

u)= E \{k'\Pd\b)\ 2 5{u-e k ). 



(i) 



(21) 



(22) 



k'>k F 



A similar result was derived and discussed earlier by Rehr 
et al 0. However, an important difference in our formu- 
lation is the presence of the projection operator P in the 
dipole matrix element. It is interesting to note that A(lo) 
contains the core electron edge singularity A(ui) « w a_1 , 
where a is the singularity index, and that fi^ ' (w) is also 
singular at the Fermi level, (i^(lv) « (uj — Ep) 13 . This 
latter singularity follows from the singular behavior of 
the overlap integral (k'\k) |l5| . We also see from Eq. j2^ ) 
that the final state rule is not strictly valid, except well 
above threshold, where P can be replaced by a (S-function 
and the matrix element reduces to (k'\d\b). We know 
from photoemission jl(| that there is little extrinsic scat- 
tering at threshold, where this limiting case should be a 
good representation of our basic approximation given by 
Eqs. (|l|) and (£f). 



Eqs. ©, (|TT]), (|T|) and @ form the basis for our anal- 
ysis of x-ray absorption. It is easy to show that Eq. ( |l2] ) 
gives a non-negative absorption cross section /i(w), as it 
should. 



B. Limiting cases 

We start our analysis of the theoretical model devel- 
oped above by discussing the two limiting cases: 1) when 
there is no extrinsic scattering (V pv = 0), and 2) when 
the core hole potential is neglected (H' v0 = H v0 ). 



No extrinsic scattering 



In this case V r 



p V 



0, 



and we can put in a complete 
set of eigenstates \<fr' n ) to H' v0 with eigenvalues E' n and 
obtain (taking 7 as an infinitesimal) 



ffcffM = E 



\(*o\K)f 



LU — LU„ 



h' 



(19) 



2. No core hole 

When the core-hole potential is neglected, H' vQ = H v o, 
h' = h, and we have (taking 7 as an infinitesimal) 



1 



uu - (H v0 - E Q ) - h - Vp V + 27 



$0 



(23) 



Now g c fi(uj) is equivalent to a standard Green's function 
g{tu)- The k\ki representation of g(u)) is 

(fci|5(w)|Afc) = 
1 



'fei 



u - (H v0 - Eq) - h - Vp V + i"f 



$0} |0) , 
(24) 

where |$o) |0) is an eigenfunction of the full Hamiltonian 
H v o + h + Vpv since V pv |0) = 0. We can express g in 
terms of a spectral function A(lu) 1 



Ei 



' (fci \A(uj')\k 2 )duj' 

LU — LU 1 + 27 

uj — h — 1j(lu) 



(25) 
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This limiting case gives a theory very similar to one- 
electron theory, but with an additional complex, energy 
dependent one-electron potential S(w). If we approxi- 
mate by a constant — iT, which is equivalent to a 
Lorentzian line shape for A(ui), we recover the conven- 
tional XAFS result, in which extrinsic losses are repre- 
sented by a mean free path Xk ~ k/T term, i.e., with 
a factor exp(— R/Xk) in each multiple-scattering path of 
length R. In general £ (u>) has structure at energies away 
from the quasi-particle energy, giving rise to satellite ef- 
fects. 



C. Quasi-boson representation 

We now turn to the general case, in which all three 
potentials V pv , V pc and V vc that couple the three sub- 
systems - the photoelectron, valence electrons, and core 
electron - are nonzero. To handle this we introduce a 
quasi-boson model Hamiltonian, 



Hvo 

V vc 



E^„4a„, 



ti 



k>k F 



v*« = £ [n r uX + (n n lfe2 )xRc fc2 . 



(26) 
(27) 
(28) 



This model together with Eqs. (|11|), (|12() and (gj) define 
the set of approximations that we use in this work. The 
potential V pc never appears explicitly since h\&o) = 0, 
and thus we do not have to worry about the transform 
between the h and h! states. 

The quasi-boson model has been discussed e.g., in Ref. 
H and [I?]. The essence is that the electron-hole type 
excitations are represented by bosons a n with energies 
Lu n , and the electron-charge fluctuation coupling is rep- 
resented by a term linear in the boson operators, as in 
Eq. (^8|). This is analogous to the usual electron-phonon 
coupling. Eq. (^) is a special case of Eq. (p8| ) with c' k c k2 
replaced by bb^ and a minus sign, because the core hole 
potential is attractive. The quantities V n are fluctua- 
tion potentials corresponding to excited states n. The 
V n can be obtained e.g., from an RPA type dielectric 
function p6| . With this simple model Hamiltonian we 
can solve explicitly for the relation between the ground 
states of H v q and H' v0 , i.e., 



3- 5 K), s 



T/n 



£ 



v, 



1,1, 



where aj, belongs to H' v0 = 5Z„w£a„a n - Expanding to 
second order in the coupling functions V n , we obtain 



n ^ J 



where (cf. Eq. (pi)), 



w-iH^-EQ-h'-V^ + i-y 
1 



E (lu) + ij ' 



(30) 



is the damped Green's function calculated in the presence 
of a core hole potential. With the above result for g e ff(u>) 
we have achieved our goal of expressing fi (w) in Eq. ( p^ ) 
as an expansion in one-particle Green's functions, thus 
avoiding the calculation of correlated many body final 
states. In the next section we take the further step of 
making a MS expansion of the Green's functions. We 
note that the limiting expressions in Sec. II. B. all come 
out the same if we had chosen to start with the quasi- 
boson model. 



D. Qualitative discussion of 2nd order expression 



While our basic approximation in Eqs. (|12|) and fllS| ) 
gives positive absorption fi(u>), there is no guarantee that 
the individual terms in an expansion in powers of the 
coupling functions V n should be positive. The coupling 
strength may be gauged by the value of the dimensionless 
coefficient a defined above. For electron gas models of 
solids, counting only plasmon modes, its value is typically 
is 0.2—0.4, and hence is quite strong. For values of k close 
to the Fermi surface the extrinsic quasiparticle strength 
Z is about exp(— a) » 0.8 — 0.7. Such values for Z are 
typical for most solids. The fact that effects of order 
a are not generally observed in XAS can be viewed as 
empirical evidence for strong cancellation effects among 
the various losses. However, given these large values of a, 
it is not surprising to encounter some non-physical effects 
in numerical calculations, such as small regions where the 
spectral function can become negative, which are due to 
the neglect of terms higher than 2nd order in the theory. 

Summarizing our second order expression, and chang- 
ing the definition of g e s in Eq. ( p"8| ) by taking out the e~ a 
factor, we have for the absorption spectrum, 



-St 



(^o - E'o) -h'-V f 



pi, 



(i {uj) = Im (b\dP ge e(w + E c )Pd\b), (31) 

7T 
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where 

9ett{w) = g qp (uj) + g ex tr(u) + gintr(u) + g m ter(u), 



+ g ex tr(u) = g(u), 



gintr{&) 



— 9(w-«n), 



V, 



9inte 



H = -2£^ fl (o;-w n )V" fl (u;),(32) 



write h' = t + Vq + V sca tt, where t is the kinetic energy 
operator, V ' the potential in the central cell with its core 
hole, and V sca tt the total scattering potential from all the 
neighboring cells (excluding inelastic losses). In the MS 
expansion we can e.g., use ho = t + V ' — iT(uj) — ij to 
obtain the propagator for the central absorber and then 
treat V sca tt as the perturbation. 



From Eqs. (J31|) and (|32[), the quasi-particle contribu- 
tion to the x-ray absorption is 



account respectively for the quasiparticle term, the ex- 
trinsic and intrinsic loss satellites, and the interference 
between them. To handle the one-particle propagators 
g(u>) we assume that T,(u>) is diagonal in a representation 
with eigenfunctions \k') of h' which according to calcula- 
tions with the GW approximation for the self energy, is 
not too bad fl7j| . For simplicity we now drop the prime 
on k' and write 



(k\g(u;)\k)=g(k,uj) = 



1 



UJ — £fc — E(fc, lo) + i'f 



(33) 



With k fixed (k \g(uj)\ k) as a function of uj has a quasi- 
particle peak and some more or less pronounced satellite 
structure. For uj near the quasiparticle peak we obtain 
an asymmetric lineshape, 



](k,uj) « (k \g qp {u))\ k) 



uj - Ek + iTi 



(34) 



where Ek is defined from Ek = fc 2 /2 + ReS (k, Ek) , Tk = 
|ImS(fe,£; fc )| and Z k =[l- dZ(k,Lu)/duj}~l Ek . We now 
make an on-shell approximation, defining functions of uj, 

Z(uj) = Z k , AE(uj) = AE k , r(w) = r fe , 

with the relation between k and lo given implicitly 
through 



u = E k . 



(35) 



Strictly speaking E is defined from h + S = t + Vu + S, 
where t is the kinetic energy and Vr the Hartree poten- 
tial. Often we would like to use e.g. Vlda rather than Vh 
to generate basis functions. We then have to replace E 
by E — Vlda + Vh in our expressions. With this on-shell 
approximation we have 



(k\g qp (cu)\k) « 



Z(oj) 



cj-h' - AE(lj) + iT(uj) + 27 



(36) 



and since the \k) are eigenfunctions of h! , the operator 
g qp {oj) becomes 



Z(oj) 



uj-h' - AE{lo) + iT(u) 



(37) 



With this form for the quasi-particle propagator g qp (uj) 
we can use the MS expansion to treat the XAFS p0| . We 



M 9P (w) 



-Im (b\SPg qp (u)Pd\b), 



(38) 



where the propagator g qp (uj) is given by Eq. <pl\). This 
is similar to the standard one-electron formula for the x- 
ray absorption fj,^> (uj) with mean free path effects from 
the damping parameter T(uj), except for a complex am- 
plitude factor Z(lo) — exp(i<fi)\Z(uj)\ with a many-body 
phase shift <f) and a wave-function overlap reduction by 
the factor e~ a . With T and a both zero and Z = 1, 
the quasi-particle XAS /i qp (uj) becomes identical to the 
one-particle absorption //^(w) in Eq. (g2|). 

To evaluate the total absorption fi(u>) including intrin- 
sic losses and interference, we can, correct to terms of 
2nd order in the fluctuation potentials V n , replace g(uj) 
by g qp (uj) in gi n t r (uj) and gi n ter (<*>). For definiteness we 
define the extrinsic satellite function g ex tr(i^) as the dif- 
ference between the full propagator g (uj) and the quasi- 
particle propagator g qp (uj) (see Appendix A). As noted 
above, the projection operator P in Eq. ( |3l| ) and Eq. ( g^ ) 
is necessary to account for edge singularity effects, but 
does not significantly affect the fine structure. In Ap- 
pendix A we also show that the many-body expression 
for the x-ray absorption /i(uj) can be expressed as a con- 
volution of an effective spectral function A c g (uj, uj'), and 
the quasi-particle absorption \i qp from Eq. (|38|), i.e., 



H{u> - E c ) = I dui'A cS (uj,uj') fi qp (u> - uj'), 



(39) 



where uj' is the excitation energy. 

To evaluate our theory numerically in real systems is 
a heavy undertaking, and we will here only carry out 
some rough estimates which illustrate the theory and 
yield non-negligible corrections to the usual XAFS pro- 
cedure [0. To simplify these calculations we rely on 
electron gas theory within the plasmon-pole approxima- 
tion to evaluate the various contributions to the effective 
spectral function, and then use these results to estimate 
the corresponding contributions to /Lt(w). Since our aim 
here is only to carry out a pilot study, which is at best 
semi-quantitative, such an approximate model seems ap- 
propriate. However, our approach is more general, and 
could be refined at the expense of much heavier calcula- 
tions. 
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III. MODEL CALCULATIONS 

In this section we present electron gas model calcu- 
lations based on the plasmon-pole approximation (Ap- 
pendix B) for the various contributions to the effective 
spectral function A e g. They can be represented as a 
sum of quasiparticle, interference, intrinsic, and extrinsic 
satellite terms defined in Appendix A, i.e., 

A cS {lo,lo') = [1 + 2a{Lu)]5(u') + A sat {Lu,Lu'), (40) 

where 

A sat (to,u') = A extr {uj,uo')+ A lntr {uo,uo') 



■24t(wV)- 



(41) 



Since we make comparisons to XAFS experiments for fee 
Cu metal, we have set r s = 1.80, which corresponds to 
the mean interstitial electron density. For this density, 
the dimensionless constant a relating the strength of the 
electron-plasmon coupling to the plasmon excitation en- 
ergy is 0.31. Near threshold the net weight of each of 
the contributions A extr , Af^| er and Ai ntr is equal to a 
and their shapes are similar, so that the sum of all of 
these contributions tend to cancel. Also near threshold, 
Z exp(— a) and the interference contribution to the 
quasi-particle peak a(u) a. Thus the net strength of 
the main peak at threshold in our 2nd order theory is 
Zexp(-a)(l + 2a) -> 1 + 0(a 2 ). 

As noted above, the asymmetric quasiparticle spec- 



trum A qp (k,Lu) = (— 1/tt) lmg qp (k, ui) 



Eq. 



gives rise to a net reduction in the XAFS amplitude as 
well as an additional phase shift compared to one-electron 
theory. These effects are due to the behavior of the com- 
plex renormalization constant Z(u>), which gives A(uS) 
an asymmetric Fano lineshape. In Fig. [jl we plot the 



modified quasi-particle spectrum A™° (k,uJ), where long 
range contribution from the imaginary part of Z^ is cut 
l™ od , (cf. Appendix A and Eq. jAl8|) ) i.e., 



A™ d (k,u + E k ) = - 



1 T k ReZ k - we-( u / 2u ») 2 ImZi 



(42) 



For comparison we also show the total extrinsic spectral 
function including both quasiparticle and satellite parts. 
Note that A™° d (k, u>) has nothing to do with the different 
contributions to A c r. The real and imaginary parts of the 
renormalization constant Z(uS) are plotted in Fig. |[ 

The intrinsic contribution Ai n t r {uj,Lo') to A e ff(u,uj') is 
independent of lo and gives a well defined satellite struc- 
ture peaking at an energy u) p away from the quasiparticle 
(Fig. ||). Although Ai n t r turns on sharply, this singular 
structure is suppressed by broadening and interference 
terms as described below. The interference between ex- 
trinsic and intrinsic losses results in a net shift of spectral 
weight away from the satellite and to the quasiparticle 
peak, overall spectral weight being conserved. The rough 
cancellation of the satellite terms due to interference is 
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FIG. 1: Quasiparticle spectral function A^° (k,io + Ek) in 
Eq. (jl2|) plotted vs u for different values of the wave number k 
(solid line); and the spectral function (cf. Eq. ( |3^ ) (— l/7r) Im g 
(dashed line). All parameters are calculated using a plasmon 
pole dielectric function for an electron gas at the mean inter- 
stitial electron density in Cu r s — 1.80. 




20.0 



IE, 



FIG. 2: Real (solid line) and imaginary (dashed) parts of 
the quasiparticle renormalization constant Z(oj) = 1/(1 — 
d'S/dcu), calculated at the quasiparticle peak, for the GW 
plasmon pole self-energy E of an electron gas at r s = 1.80. 
The sharp structure occurs at the onset of plasmon excita- 
tions. 



clearly illustrated by the behavior of A\^ ter in Fig. ||. 
Note that the interference satellite amplitude is maxi- 
mal near threshold and slowly decreases with increasing 
energy over a range of several ui p . 

The behavior of the extrinsic satellite spectral func- 
tion is illustrated in Fig. ||. In the plasmon-pole, electron 
gas model used here, the extrinsic satellite spectral func- 
tion sometimes exhibits a complicated structure. Close 
to the Fermi energy, the structure simplifies and con- 
sists of a peak at an energy about oj p above the quasi- 
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FIG. 3: Intrinsic (solid line) and interference satellite 
(dashed) spectral functions Ai ntr (uJ, to 1 ) and A°%l er (uj,uj'), 
plotted against uj' for selected values of uj. Note that the 
interference contribution is negative and tends to cancel the 
intrinsic satellite. These quantities are both calculated us- 
ing the plasmon-pole dielectric function for an electron gas at 
density r s — 1.80, and broadened by a Lorentzian of width 
0.2uj„. 



FIG. 4: Extrinsic satellite spectral function A ex tr(uJ, uj') ob- 
tained as described in the Appendix for selected values of uj, 
vs uj' . The curve for uj — 0.7J5f and uj = 1.5Ef lie before 
the onset of plasmon excitations. The result for uj — 1.5Ef 
illustrates the anomalous structure seen in this region, and 
that for uj = 6Ef lies well beyond the onset. All results are 
based on a plasmon-pole dielectric function with r s = 1.80. 



particle energy uj and a smooth structure that falls off 
gradually with increasing energy. For u> near the on- 
set of plasmon losses, there is still a pronounced satel- 
lite peak, but there is also an additional "anomalous" 
structure near the quasiparticle peak. Indeed, it seems 
ambiguous whether the structure close to the quasipar- 
ticle energy should be considered as part of the satellite 
or the main peak, as the structure accounts for a sub- 
stantial portion of the extrinsic weight [1 — Z(uj)] that 
is not included in the quasiparticle peak. This indicates 
that the anomalously low and singular behavior of Z(uj) 
in this region is partly due to the singular structure of 
the plasmon-pole approximation and largely an artifact 
of the ad hoc method used to separate the main peak and 
the "satellite" spectral function. Above the onset of plas- 
mon losses, the anomalous structure disappears, and is 
replaced by a small tail which extends to the vicinity of 
the quasiparticle peak. Moreover, as the quasi-particle 
energy increases, the extrinsic satellite weight becomes 
progressively smaller. 



IV. IMPLICATIONS FOR XAFS 

We can now obtain rough estimates for the effect of ex- 
trinsic and intrinsic losses and interference on the XAFS 
spectrum. In the usual MS theory Jl(J], the XAFS spec- 
trum (w) is a rapidly varying energy dependent factor 
in the one-particle expression for the x-ray absorption, 



M W(a;)=^H[l + X W( W 



Mi 



(43) 



where /Iq 1 ' is the generally smooth absorption from the 
central atom alone, in the absence of MS. The conven- 
tional (broadened) one-particle absorption /Z 1 ) (uj) is ob- 
tained from Eq. ("J)), with g e s(u) = 1/[uj — h — iT(uj)], 
i.e., with a damped one-particle propagator with mean 
free path effects taken into account in terms of T(ui). 
This one-particle MS theory is generally in good agree- 
ment with experiment. However, there remains a resid- 
ual discrepancy of about 10% in overall XAFS ampli- 
tudes and a systematic shift in peak positions compared 
to experiment. This shift is only partly accounted for by 
including the real part of an electron gas self-energy in 
the one-particle propagator. 

In the present theory, the many body effects of losses 
and interference can be represented as in Eq. (^9|), i.e., as 
a convolution of fj, qp (tu) with the effective spectral func- 
tion A e g(u), u)') including contributions from both pri- 
mary and satellite channels, as discussed in Appendix 
A. As noted above, the difference between the behav- 
ior of fj.^ (uj) and [i qv (uj) is qualitatively minor and re- 
sults primarily from a constant wavefunction overlap fac- 
tor exp(— a), a complex renormalization factor Z(uj) = 
\Z(iS)\ exp(i(f>), and an energy shift AE (uj). The phase- 
factor (f> can be absorbed by adding <fi/2 to the central 
atom phase factor in the XAFS formula. 

To extract the XAFS x( w ) corrected for many-body ef- 
fects, we first use a similar MS factorization of fi qp (uj) = 
A*gp[l + Xgp( w )] to split off a central cell contribution 
(uj), which yields a many-body expression for the 
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atomic background absorption 

MoM = J ' duj'A cff (uj,uj')fi qp (Ld - lu') 

dw'A eS {u,uj'). (44) 

-oo 

Here yP qp {u) ~ | |P<i|6)| 2 is the absorption from the 
central atom alone, in the absence of other scatterers, 
and we have neglected the variation of ^ p (u>) over the 
dominant integration range of a/. This is usually a good 
approximation, since the one-particle atomic background 
is usually a smooth, monotonically decreasing function of 
energy. Thus one expects that the satellite structure in 
the spectral function due to many-body excitations will 
generally have minor effects on fj,o(u>), which is consistent 
with experimental observations. However the existence 
of sharp atomic resonances in the one-particle absorp- 
tion n° qp may lead to exceptions. The many body XAFS 
function x( w ) — (/•* — A*o)/mo then becomes 

xM ~ J duj'A cS {uj, uj')x qp (uj ~ u'), (45) 

where the spectral function A c g is now normalized to 
unity, 

i cff («, «/) = A cS (lo, u/)/N(u), (46) 

and7V(cj) = J du)' A eS (to,uj'). With this normalization in 
x(^), the wave function renormalization factor exp(— a) 
and the magnitude of the quasi-particle strength \Z\ can- 
cel out, while the phase <j> and the energy shift AE(lu) 
both remain. 

The importance of multi-electron excitations can be 
gauged by the net spectral weight in main peak and in 
the satellite structure from all losses. A plot of the nor- 
malized integrated satellite spectral weight, 

a sat (cj) = J du'A sat (u),u') (47) 

where A sat (u>, uj') — A sat /N(uj) is given in Fig. ||. Also 
plotted is the total weight of the primary peak, 

Oo(w) = [l + 2o(w)]/JV(a;). (48) 

Note, in particular, the slow trend of the satellite weight 
towards the sudden approximation limit [aexp(— a)] with 
increasing energy over a range of a few hundred eV. In 
these plots we have lumped the contributions to the spec- 
tral function that lie below the plasmon onset into the 
primary peak. The anomalous behavior of the weights 
near the plasmon onset energy is due to the ambiguity of 
separating the satellite and quasiparticle contributions, 
and does not lead to singular structure in the overall ab- 
sorption. 

The net effect of the convolution over a normalized, 
positive spectral amplitude A e g(u;,u/) of Eq. ( p5| ) is 




u/E F 

FIG. 5: Normalized total spectral weights of the net pri- 
mary peak (solid) and satellite terms (dashed) including in- 
terference, i.e., the integrated quantities a aat {iu) and clq(uj) of 
Eq. (^) and (^). The anomalous structure at low energies is 
an artifact of the method used to separate primary and satel- 
lite terms in the spectral function. All results are obtained for 
the model of this paper with r s — 1.80. The horizontal line 
represents the high energy limit of the satellite weight (i.e., 
the sudden approximation). 

clearly a decreased XAFS amplitude and a phase shifted 
oscillatory signal compared to the one-particle XAFS 
■ In the single scattering approximation the os- 
cillatory energy dependence of Xgp(^) enters primarily 
through the complex exponential exp[i2k(u))R], where R 
is an interatomic distance and k(u) — ^/2uj is the pho- 
toelectron wave vector. The result of the convolution 
can be written in terms of a complex amplitude factor 
Sq(u,R) = \Sq(uj, R)\ exp(?$(w, R)), which is given by 
an energy dependent "phasor sum" over the effective nor- 
malized spectral function, 

S%(u,R) = [ dLu'A cS (uj,iu')e l2 ^-"">- k ^ R . (49) 

The qualitative behavior of Sq(lu,R) can be understood 
as follows: At very low energies compared with the exci- 
tation energy u> p , the satellite terms strongly cancel so 
A(cj,u') ps S(oj - u/) and hence, Sq(uj,R) -> 1. At 
high energies, the sudden approximation prevails, and 
A ps A qp + Ai n tr, which has a strong satellite structure. 
However, the phase difference 2[k(u> — u>') — k(u>)} be- 
tween the primary channel and satellite becomes small 
at high energies (u/ ^> uj p ) and hence also Sq(u},R) — > 1. 
At intermediate energies, however, the value of Sq(oj,R) 
has a minimum. A plot of the magnitude and phase of 
Sq(lj,R) for our electron gas model is given in Fig. j|, 
for the first neighbor distance of Cu metal R — 2.55 A. 

In order to compare these results with experiment, we 
isolate the first shell of the experimental EXAFS sig- 
nal x( w ;-R) by Fourier filtering over the range 1.75 < 
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5.0 10.0 15.0 20.0 

EXAFS Wavenumber ^(A" 1 ) 

FIG. 6: Upper curves: magnitude \Sq\ (solid) of the XAFS 
many-body amplitude reduction factor due to all inelastic 
losses, as calculated using a phasor summation over the total 
spectral function for the model in this work, and plotted vs 
the EXAFS wave number k = y2(u — Ef) for the first neigh- 
bor distance R = 2.54 A in Cu. Shown for comparison are 
the experimental (dots) and theoretical (dashes) amplitude 
reduction obtained by Fourier filtering (see text). Note that 
the largest discrepancy between theory and experiment occurs 
where |x (k, R)\, the Fourier filtered and back-transformed 
complex XAFS amplitude, is small and experimental noise 
dominates. Lower curves: phase $ of So(k,R) (solid), and 
comparison, the many body phase shift of the asymmetric 
quasiparticle peak (dashes) i.e., (j> = tan _1 [ImZ/ReZ]. The 
sharp structure near k ~ 2 is an artifact of the sudden plas- 
mon onset in the plasmon pole model used here. 



R < 2.80 A in position-space R (conjugate to 2k), with 
a smooth sine window, and then back-transforming to 
fc-space. For this comparison, we use the usual EXA FS 
convention for the wave number k = \J2(E — Ep), as 
measured from the threshold Fermi energy. Our estimate 
of the experimental Sq(u>,R) is then given by the ratio 
of this back-transformed experimental first shell XAFS 
signal x(w, R) to a similarly Fourier filtered and back- 
transformed theoretical first shell signal x^(uj,R). The 
latter is obtained from ab initio XAFS calculations using 
the FEFF8 code ||. The FEFF8 calculations include 
only extrinsic losses, i.e., the mean free path loss calcu- 
lated from a Hedin-Lundqvist plasmon-pole self-energy 
model. The results, from both theory and experiment 
are also plotted in Fig. |[ Due to the Fourier filtering, 
fine details of the theoretical phasor sum for S^(uj,R) 



in Fig. are lost. Also plotted, in Fig. |. is the back- 
transformed XAFS amplitude |x(u;(fc), R) for the first 
nearest neighbor of Cu. This illustrates both where the 
amplitude reduction is important in analysis, and also 
where x is small, and hence the Sq(oj,R) extracted from 
experiment may have significant errors due to experimen- 
tal noise. Given the rough, electron-gas approximations 
used in our model calculations, the overall agreement 
with experiment for Cu metal is reasonably good. This 
result also suggests that the conventional procedure of 
approximating Sq(uj,R) by a constant Sq « 0.9 is not 
unreasonable. The biggest discrepancies are at low ener- 
gies, and are likely due both to experimental noise and 
to the approximation used for the mean-free path at low 
energies in FEFF8, which often has too much loss. We 
have also plotted the many-body correction to the XAFS 
phase, which varies by about +0.2 radians over the XAFS 
experimental range 3 < k < 20 A -1 . By comparing with 
the phase of the renormalization constant Z (Fig. |2|), one 
sees that much of the phase shift arises from the asym- 
metry of the quasiparticle peak. The sign of the phase 
shift is consistent with a reduction in the strength of the 
self-energy due to cancellation effects. The approximate 
linear variation of the phase with k can lead to errors 
in distance determinations from XAFS measurements of 
about SR = A$/2Afc ps +0.006 A which is comparable 
to systematic errors typically encountered in experimen- 
tal XAFS analysis. 

Finally we plot in Fig. [?] a comparison with experiment 
of the full absorption result /i(w) obtained by convoluting 
the spectral function A e g with the quasi-particle result 
fi qp {u>). These calculations were carried out using full- 
multiple scattering calculations for a cluster of 300 Cu 
atoms using the FEFF8 code. To compensate for errors 
in our 2nd order expression for the spectral function, A c g 
was normalized to unity in this convolution; also since 
fiqp is real, the asymmetry in the quasiparticle peak in 
Eq. ( p^ ) was represented by a term PlmZ/u; (where P 
denotes the principal part) rather than a complex phase 
factor expi(f>. Although it is not clear how reliable such 
a calculation is at the edge, given the simplicity of our 
model calculations, the agreement with experiment does 
markedly improve. 



V. CONCLUDING REMARKS 

We have developed a semi-quantitative theory for the 
effects of extrinsic and intrinsic losses and the interfer- 
ence between them on x-ray absorption spectra. The 
theory is based on a quasi-boson representation for the 
excitations, and perturbation theory to 2nd order in the 
electron-boson coupling. These losses lead to an asym- 
metric broadening of the main quasi-particle peak, plus a 
broad energy-dependent satellite in the spectral function 
due to the boson excitation. We find that the interference 
terms strongly suppress the satellite terms and enhance 
the main quasi-particle peak. These results thereby ex- 
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FIG. 7: Comparison of the full calculated many-body XAS 
/i(a>) (solid) obtained by convoluting the spectral function A e ff 
for the model of this paper with the one-particle XAS /i' 1 ' (o>) 
(dashes) calculated using the FEFF8 code, and the result from 
Cu experiment (dots). 



plain the general weakness of multi-electron excitations 
in XAS. We have applied the theory within the electron 
gas approximation to estimate the many-body correc- 
tions to the XAFS. By using a phasor summation over 
the spectral function, the theory yields an approxima- 
tion to the reduction in the XAFS amplitude beyond the 
usual extrinsic, mean-free-path, together with an addi- 
tional many-body phase shift. Pilot model calculations 
based on the electron gas approximation and our for- 
malism are in semi-quantitative agreement with XAFS 
experiment. 

Our results illustrate a striking difference with those 
for losses in photoemission. In core electron photoemis- 
sion for metals there is also a strong extrinsic contribu- 
tion to the losses for all energies, and a strong interfer- 
ence effect between extrinsic and intrinsic losses up to 
energies in the keV region pq |. In contrast, we found the 
extrinsic losses to have a very small influence on XAS. 
This at first seems to be a paradox, since photo absorp- 
tion is usually considered to be proportional to the pho- 
toyield. This is indeed true in one-electron theory, and 
is also observed to be approximately true experimentally 
in many cases. There is, however, no theoretical justi- 
fication for assuming this to be the case in general. If 
we examine the problem formally, the XAS is related to 



the dielectric response function (i.e., the density-density 
correlation function), while the PES is given by a three 
current correlation function, and there are no simple con- 
nections between these quantities. More physically we 
can point at the fact that in PES we measure a sharp, 
well-defined photo electron; i.e., we can describe PES by 
the golden rule with well defined final states. In XAS the 
photo electrons are never measured, and photo-electron 
states are not sharp, they are decaying states with a fi- 
nite lifetime. Further all electrons which are photo ex- 
cited do not leave the solid, and some PES current which 
leaves the solid is due to secondary electrons which can- 
not come back and interfere at the photo excited atom. 
We never consider these secondary electrons since we do 
not allow our quasi-bosons to decay. If we could put 
in well-defined intermediate states with a photo electron 
in the density-density correlation function we would in- 
deed have the same matrix elements that enter in PES, 
but that is not possible except perhaps approximately, 
very close to threshold when the quasi-particle life-time 
is long. 

Finally we would like to comment on similarities and 
differences between our approach and a model studied 
by Schrieffer fljfl . That work considered an Anderson- 
Newns like model for PES from a valence level on an 
adsorbed atom, coupled to surface plasmons. Schrieffer 
studied photoabsorption, and by cutting the polarization 
diagram, he identified a PES final state and obtained a 
perfect square expression for the PES current, including 
interference effects. His results agree precisely with what 
we have in our model both for XAS and when applied to 
PES. However the results for XAS cannot be expressed 
as a perfect square, as needed to have a close correspon- 
dence with PES. This is due to the differences in the 
signs of the imaginary parts, as can be seen by comparing 
Eq. (7), (13) and (14a) in Ref. |l9[ Thus photobsorption 
and photoyield arc only approximately related, even for 
a finite system where the quasi-particle aspect does not 
enter. 
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APPENDIX A: DERIVATION OF THE 
EFFECTIVE SPECTRAL FUNCTION 

In this Appendix we derive expressions for the differ- 
ent parts in the spectral weight function A e R = A qp + 
Aextr + Aintr + Ai nter that appears in the convolution 
expression for the x-ray absorption in Eq. ( |39| ) . Starting 
from Eq. (31), we have 



ti(Lu-E c ) = -—y^\(b\JP\k)\ 2 Img cS (k,Lu), (Al) 



provided g e e (see Eq. p2| ) is diagonal in the one-particle 
index k. From GW calculations we know that g and 
g qp are approximately diagonal and this is hence 
a reasonable approximation for all contributions to g e g 
except 3i„ter: which requires a special treatment as de- 
scribed later in this Appendix. We first introduce the 
real spectral weight function A c ff(k,uj) by representing 
Im<7eff(fc,w) as 

Img cS (k,Lj) — J ' duj 1 'A cS (k,uj')Img qp (k,Lj - lj'). (A2) 

Next we make the on-shell approximation of replacing k 
in Aes(k, u>') by k = k{oj — o>'), with k(tv) defined in Eq. 
(||). With A eS (u>, lo') = A cS (k (lj - a/) , J) depending 
only on lo and lo' , we can perform the summation over k, 
and obtain 



(i(u> - E c ) = I dw' A e $(u,u') Hgp(u> - lo'), 



(A3) 



where fi qp (lo) is defined in Eq. p8. Thus the quasi- 
particle contribution to A c s (lo, u7) is simply a delta- 
function, A qp (uj,uj') — 5(lo'). 



For the intrinsic contribution, it is clear from Eq. (|3^ 
that A>ijit r is simply a sum of energy-shifted <5-functions 



A mtr (Lu,tu') = Y<{ — 



5{J 



(A4) 



which is independent of lo. In the electron gas, plasmon 
pole approximation the explicit expression for Ai n t r can 
easily be found (Ref. [li], pp 655), 



A intr (lJ,Lj') 



L0 2 p e(L0' - W p ) 



(A5) 



where q(co r ) is a solution to to q — to' . The intrinsic 
spectral function contributes only at energies near to p 
above the electron quasi-particle energy, and gives a well- 
defined satellite. 

The interference contribution to g c g in Eq. ( |32| ) con- 
tains a product g(ui—uj n )V n g(uj), and since V n can trans- 
fer momentum, is not diagonal in k space even when g is. 
One can evaluate ginter{u), e.g., by doing MS expansions 
of the two propa gat ors, but this does not yield a result of 
the form in Eq. (EM) . To force it into that form we make 
some further approximations. First we take the fluctua- 
tion potentials as plane waves, and in a plane wave basis 
we then have (k + q'\V q \k) = V^S^*. Also taking the 
Green's functions as diagonal in a plane wave basis, Eqs. 
(|3~i"|) and ||| give 



Winter 



2e~ 



-Im 



kq 



l^L(6|dtp|A- 

UJ a 



q) 



x(k + q\g (u - U) q ) \k + q) (k\g (w) \k)(k\Pd\b). (A6) 



Taking the quasi-particle approximation in Eq. (34) for 
g (to) and neglecting the (/-dependence in the dipole ma- 
trix element, we have 

Winter M = — ImV^- \(k\Pd\b)\ 2 



TV 



kq 



+ Ek+q — Ek 
1 



W — Wq — Ek+q + iTk+ q to - E k + iY k 



,(A7) 



where we have made the approximation that Tk+ q — Tk 
in the denominator of the prefactor before the term in 
brackets. We now have a difference between two Green's 
functions, and provided we treat the prefactor of this 
difference as a real number (e.g., neglecting the imaginary 
part of Z) we can again write Winter ( w ) as a convolution 
with /Xgp, i.e., 



Winter (w) 



a(u))ii qp (u)) 



,(A8) 
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with 



IVI 2 



LU q Ul q + Efr+q — E k 



(A9) 



k-k(uj) 



inter 



(uj,uj') = 



</|2 



\ Z k-q\ 



1 



us q + E k - E k - 



X 5(u>' — UJq) . 



k=k(uj—uj') 

(A10) 



To further simplify the calculations, we assume that 
Zk+q ~ Z k . Thus we have separated the resulting spec- 
tral function into a term proportional to a(u>) that adds 
to the quasi-particle peak, and a term proportional to 
Afnter that contributes to the satellite structure. Ap- 
proximating Ek by k 2 /2, we obtain 

a(us) = pl / dus' 



Tirk 



x log 



^'(Wp 

LUp + 2lu' + fcg 



,/\2 



up -r 2u>' — kq 
where q = v2u/ and k = \f2uu. Similarly 



(All) 



4 s 



(u>,u>') = 



\Zy\u! 



P 



x lo. 



4irk'(w' - lu p )lu' 2 
ojp + k'(us)q 



where q — y/2(u>' 



ujp — k'(us)q 
and k! = ^2(us - us') 



lu' -lu p ), (A12) 



For the extrinsic contribution, we have from Eq. ( |A2| ), 

Qextr 

(k,Lu) f A extr (k,us')dus' 



Im 



Z k 



Ek + iTk 



= 0. (A13) 



The integral gives an analytic function of lu in the half- 
plane above the line lu — Ek — iTk ■ If we also assume that 
9extr (k, lu) — g (k, Lu)—g qp (k, lu) is analytic in this region, 
we can make a shift in the origin for us by Ek — iTk + iS 
(5 — > + ) without crossing the analyticity line. This gives 



Im 



gextr (k, us + E k - iT k + iS) 



Zk 

Aextr {k,Lu') 



dLu' 



lu — us' + i5 



0, (A14) 



from which we obtain 
1. 



A extr (k,us) 



-Im 



7T 



9extr{k, lu + E k - iT k + iS) 



(A15) 



The analyticity of g extr (k,us) follows if it can be de- 
scribed by a sum of discrete poles or by a contour in- 
tegration below lu — Ek — iTk- This is a reasonable ap- 
proximation since the satellite peak is broader than the 



quasi-particle one. If for g qp (k, us) we take the expression 
valid for us close to the quasi-particle energy 



g qp (k,u>) = - 
us 

and use it for all us, we have 

1 



Ek + iTk 



(A16) 



Qextr 



(k,us) 



u-e k -E(k,us) us-E k 



iT k 



(A17) 



Since Zk has a fairly large imaginary part, this leads to a 
substantial Fano type asymmetry and substantial (can- 
celling) negative contributions in the spectral functions 
for g ex t r and g qp . The expression to use for g qp is not well- 
defined except very close to the quasi-particle energy. To 
avoid a negative spectral density for the extrinsic satel- 
lite we introduce a Gaussian cutoff to the imaginary part 
of Z k in lmg qp , i.e., 



Im g qp (k, us + E k - iT k + iS) 

-5Re Z k + we-W^'lm Z k 
us 2 + S 2 



(A18) 



If we further replace Zk by \Zk\ in Eq. (A15), we obtain 
for lu >> S 



~(k, us) 



1 



nZ k 



[T k + ImJ:(k,LU + E k 



xl 



[us + AE k - Re E(fc, us + E k )Y 



[T k + lmY,(k,uj + E k )Y 



ImZk 



(A19) 



where AEk = Ke^(k, E k ). We have neglected the iT k 
in the argument for S, and we have only considered 
us >> S since g extr (k,us) is small around lu = when 
the true complex quasi-particle energy is used. This re- 
sult can be compared with the first order result from 
Eq. (A2), i.e., the result obtained by taking Tk — > 0. 
With T k = and neglecting the phase in Z k we have 
Im g qp (k, lu) — — it \Z k \ 8{uu — E k ), and thus from Eq. ( [A2] ) 
{k, lu - E k ) = -1/ (it \ Z k \) lmg extr (k,uj). This is 
the same result as in Eq. (A19) with = 0. Since 
the GW approximation puts the extrinsic satellite at a 
slightly different position than the interference and in- 
trinsic satellite terms, (which leads to numerical prob- 
lems such as small regions where the spectral function is 
negative) we shifted A ex ± r to make the peaks coincide. 
In summary we have found that 



fi(us-E c ) = Jdus'Aes(ui,us')ii q p(u) — u)'), (A20) 
where 

A cS (lu,lu') = [1 + 2ci(us)]5(lu') + A extr (us,ui') 

+ A mtr (u>,uj')-2A% t ter (Lu,Lu'). (A21) 
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APPENDIX B: 



PLASMON-POLE, 
GAS MODEL 



ELECTRON 



In this Appendix, we briefly outline some properties 
of the electron gas model used in our calculations. In 
order to estimate the many-body effects on XAFS spec- 
tra, we need to make some simplifying approximations. 
By choosing to work with an electron gas model, many 
of the formulae in our model can be found analytically, 
greatly simplifying the calculations. We further choose 
to use a plasmon pole dielectric function. Although the 
model exhibits some non-physical singular structure and 
no loss at low energies, it never the less gives mean free 
paths and self-energy shifts in reasonable agreement with 
experiment. Under this approximation, the fluctuating 
potentials V n are plane waves 



where u>th is the threshold for plasmon excitation, and q\ 
and q2 are limiting values of the inequalities uj q + (q — 
k) 2 /2 — uj < and uj q + (q + k) 2 /2 — uj > 0, and hence 
are solutions to a cubic equation, 



kq 3 + ( uj + a - ^k 2 ) q 2 + (k 3 - 2ujk)q 



2 2 
UJ p -UJ 



;k 2 - 



4 



= 0, 



(B3) 



with the constraints, min(qi, q2)—0, and max(gi, <72)=0 
satisfies uj q = uj + k 2 /2 — k 2 F /2. The real part is then 
obtained by a Kramers-Kronig transformation, 



V n (r) = Vi(r) 



V q e iqr . 



In the case of coupling to the core hole (at r = 0) 
this yields VA — V^. According to the plasmon pole 
model §0f 



ReE(fc,w) = V ex (As) 



P f Im£(fc,u/)cL/ 



(B4) 



2ne 2 uj; 
q 2 uj q tt 



1/2 



(Bl) 



where f2 is the system volume. For the case when the 



plasmon dispersion has the form uj q 



+ aq 2 + g 4 /4, 



the imaginary part of T,Q\y(k,uj) can be obtained ana- 
lytically, i.e., 



Im S(/c, uj) 



4A- 



log 



+ 3i 



4 



e{uj 



Uth), 

(B2) 



where 



V ex (k) 



kp 
~2~ 



k 2 



4k 



In 



k + kj, 



k — kp 



(B5) 



is the energy-independent part in E, i.e., the Hartree- 
Fock exchange energy. In the calculations presented in 
this work we have chosen the plasmon dispersion as in 
Ref. ([^O)), i.e., a = 2/3, which is the same as that used 
in the extrinsic loss calculations in the FEFF8 code. 



